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ABSTRACT 


A spray atomization and combustion model is developed based on the volume-of- fluid (VOF) 
transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and 
energy conservation laws are modeled by continuum surface force mechanisms. A new solution 
method is developed such that the present VOF model can be applied for all-speed range flows. 
The objectives of the present study are: (a) to develop and verify the fractional volume-of-fluid 
(VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas- 
liquid) free surface flow problems; (b) to implement the developed unified algorithm in a general 
purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with 
droplet dynamics and finite-rate chemistry models; and (c) to demonstrate the effectiveness of the 
present approach by simulating benchmark problems of jet breakup/spray atomization and 
combustion. Multiphase fluid flows involving free surface fluids can be found in many space 
transportation and propulsion systems such as injector atomization in Space Shuttle Main Engine 
(SSME) combustors, cavitation in liquid rocket pump operations, handling of cryogenic liquids in 
a micro-gravity environment and the operation of cryogenic propellants on board spacecraft. 
Modeling these types of flows poses a significant challenge because a required boundary must be 
applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can 
range from incompressible to high-speed compressible flows. The flow-process modeling is further 
complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, 
and their interactions. 

The major contribution of the present method is to combine the novel feature of the Volume of 
Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient non- 
iterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The 
proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and 
gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified 
predictor-corrector technique to account for the non-linear linkages through the convective 
contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume 
force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the 
droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. 
Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily 
complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as 
the baseline of the current development. 

Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and 
impinging jet test cases are investigated in the present work. Preliminary data comparisons show 
good qualitative agreement between data and the present analysis. It is indicative from these 
results that the present method has great potential to become a general engineering design analysis 
and diagnostics tool for problems involving spray combustion. 
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FLUID DYNAMICS PROBLEMS WITH MULTI-PHASE INTERFACE ARE 
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Transport Equations 
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For a given solution of a field, the governing equations can be recast as: 
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Continuum Surface Force Model 
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COAXIAL INJECTOR TEST CASES 
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COAXIAL INJECTOR SIMULATION (1 atm, Steady-State Approach) 
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COAXIAL INJECTOR SIMULATION (1 atm, Transient Approach) 
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COAXIAL INJECTOR SIMULATION (1 atm, Transient Approach) 
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UNI-ELEMENT COMBUSTOR SIMULATION (30 atm, 2-Step Model) 
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UNI-ELEMENT COMBUSTOR SIMULATION (30 atm, 2-Step Model) 
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Mach Number Contours 


UNI-ELEMENT COMBUSTOR SIMULATION (30 atm, 7-Step Model) 
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IMPINGING INJECTOR TEST CASES 
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LIQUID SHEET THICKNESS COMPARISONS (0 = 30°) 
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LIQUID SHEET THICKNESS COMPARISONS (9 = 45°) 
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LIQUID SHEET THICKNESS COMPARISONS (0 = 60°) 
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